Visualization of semileptonic form factors from lattice QCD 
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Comparisons of lattice-QCD calculations of semileptonic form factors with experimental measure- 
ments often display two sets of points, one each for lattice QCD and experiment. Here we propose 
to display the output of a lattice-QCD analysis as a curve and error band. This is justified, because 
lattice-QCD results rely in part on fitting, both for the chiral extrapolation and to extend lattice- 
QCD data over the full physically allowed kinematic domain. To display an error band, correlations 
in the fit parameters must be taken into account. For the statistical error, the correlation comes 
from the fit. To illustrate how to address correlations in the systematic errors, we use the Becirevic- 
Kaidalov parametrization of the D nlf and D Klv form factors, and an analyticity-based fit 
for the B — > ttIu form factor /+. 

PACS numbers: 13.20.Fc,13.20.He,12.38.Gc 



The past several years have witnessed considerable im- 
provement in our understanding of semileptonic decays 
of D and B mesons. Measurements have advanced in 
accuracy from 6-20% on the normalization [l|J3, 0] and 
~ 10% on the shape [1] to 1% on both Mean- 
while, ab initio calculations in QCD with lattice gauge 
theory have become realistic [1, 0, H, |^ , now incorporat- 
ing the effects of sea quarks that were omitted in earlier 
work 0, [m, [H, [H, d. In this article, we discuss how 
to present both together, so that the agreement (or, in 
principle, lack thereof) is easy to assess. 

We focus on reactions mediated by electroweak vector 
currents, leading to pseudoscalar mesons, tt or K , in the 
final state. At the quark level, a heavy quark h decays 
into a daughter quark d (not necessarily the down quark) , 
with a spectator antiquark q. Writing the decay H — > 
Plv, the form factors are defined by 

(PIV^IH) = f+{q^){pH +PP - AT + Mq^)A^^ (1) 



denote the 3-momentum and energy of the final-state me- 
son in the rest frame of the initial state. The energy E 
is related to q"^ via 



2mHE. 



(3) 



Neglecting the lepton mass, < g'^ < qf^ax — {rriH—mp)'^ 
is kinematically allowed in the semileptonic decay. 

The form factors f+{q^) and fo{<l^) are related to /|| {E) 
^nd f^{E) by 



f+{q') = {2mH)-'/' [fi\iE) + {mn ~ E)f^iE)] , (4) 



[{mH^E)f\\iE)-plME)] , (5) 



imply 



where q = pn — pp is the 4-momentum of the lepton 
system, and = {pH+Pp)-qq^ / (f' = {'m% — mp)q'^/q'^. 
Equation ^ is general, applying to K ^ -kIv as well as to 
D and B decays. For lattice QCD, it is more convenient 
to express the transition matrix element as 



{P\V^\H) = ,/^[v^^f\\{E)+py^{E)\ , (2) 
where v — pn/mH, and p± = pp — Ev and E — v ■ pp 



with Eq. 13]) understood. Equations ([4]) and 
/_^_(0) = /o(0), as required in Eq. ^. 

Two aspects of lattice-QCD calculations are important 
here. First (as in all lattice-QCD calculations), it is com- 
putationally demanding to have a spectator quark with 
mass as small as those of the up and down quarks; for 
P — TT the same applies to the daughter quark. In re- 
cent unquenched calculations, the mass of the qq pseu- 
doscalar Pqq lies in the range 0.1m|- < 'rnp_^ < mj^. Sec- 
ond (of special importance in semileptonic decays), the 
calculations take place in a finite spatial volume, so the 
3-momentum takes discrete values. In typical cases the 
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box-size L Ki 2.5 fm, so the smallest nonzero momentum 
P(i,o,o) = 27r(l,0,0)/L satisfies |P(i,o,o)l ~ 500 MeV. 

After generating numerical data at several values of 
{E,mp_J, the next step for lattice-QCD calculations is 
to carry out a chiral extrapolation, ?Tip_^ — > m^, of the 
data for and /|| [H, The chiral extrapolation 

must refle ct the fac t that the form factors are analytic 
in _E = \/p'^ + ^p, not p Note also that f± and 

/-!_ can be computed only with p =/= 0, hence E > mp 
or, equivalently, < (/nj^x- The statistical and dis- 
cretization uncertainties in [q^ , TOp_ ) start out small- 
est at q^^ Q corresponding to P(i,o,o)- A sensible chiral 

extrapolation will propagate this feature to /+(g^,m^). 
Similarly, the statistical and discretization uncertainties 
in /o(g^,mJ) are smallest near g^^x- 

When \p\a becomes too large, discretization effects 
grow out of control. Therefore, the kinematic domain 
of lattice-QCD calculations is limited to, these days, 
IpI ^ 1 GeV, with a corresponding upper limit on E 
and lower limit on q^ . To extend the form factor over 
the full physical kinematic domain, a parametrization of 
the q^ dependence is needed. 

One choice is the Becirevic-Kaidalov (BK) ansatz [Ig] 



F 



F 



aqp) 



(6) 
(7) 



where = q^/m^, {H* is the vector meson of flavor hd), 
and F, a, and (3 are free parameters to be fitted. A key 
feature of Eq. © is the built-in pole at q"^ = m%, , or 
E = —{mjjt. — mfj — rn?p)/2mH < 0, an indisputable fea- 
ture of the physical /+. Further singularities at higher 
negative energy are modeled by the BK parameters a 
and (3. A similar possibility is the Ball-Zwicky (BZ) 
ansatz IToj . which has one more parameter for /+ 
than BK. A shortcoming of these parametrizations is that 
comparisons of lattice-QCD and experimental slope pa- 
rameters can be misleading 0,[2l|, because lattice-QCD 
slopes are determined near q^ = g^^x; whereas experi- 
mental slopes are determined near q^ — 0. 

Another approach based on analyticity and unitarity 
is to write the form factors as 



1 



1 ^ 



N 



(8) 



(9) 



where (/>+,o are arbitrary, but suitable, functions, and the 
series coefficients are fit parameters. The variable 



(10) 



^\-q^/t+ + ^l-to/t+' 

where i+ = {mn + mp)^ and to can be chosen to make 
1 2: 1 small for all kinematically allowed q^. Like BK and 



BZ, Eq. dH]) builds the H* pole into /+, but this approach 
is model independent because unitarity [13, [H, [23l] and 
heavy-quark physics [2l| impose bounds on X^fel^fcP' 
J2k l^fcP; a-iid because kinematics set \z\ < 1. Conse- 
quently, the series can be truncated safely, once addi- 
tional terms are negligible compared to other uncertain- 
ties in the analysis. 

In all approaches the output of an analysis of lattice- 
QCD form factors is a fit, usually a two-stage fit of chiral 
extrapolation followed by q'^ parametrization. Clearly, 
the final fit describes a curve, and the error matrix of 
the fit parameters describes an error band. Nevertheless, 
lattice-QCD results usually have been plotted as a set 
of points with error bars at fiducial values of q^ (or E). 
These points evoke the underlying discrete nature of the 
3-momentum p but, in general, the chosen values of q^ (or 
E) have nothing to do with the original discrete values 
of p. A plot with a curve plus error band exhibits the 
same information, while giving a visually superior sense 
of the correlations between points on the curve. 

The experimental measurements of /_|_(g^) come from 
counting events in bins of q^ and removing coupling and 
kinematic factors. The analysis inevitably entails some 
fitting, to correct for acceptance, etc., but the postfit bins 
of q^ faithfully mirror the input to such fits. 

If one would like to compare the calculations with the 
measurements, it is appealing to represent one as a curve 
with error band, and the other as points with error bars. 
Bearing the foregoing remarks in mind, it seems natural 
to draw the curve for lattice-QCD calculations. A few 
years ago, we prepared illustrative plots for D Klv 
with the Fermilab-MILC m, B lattice-QCD calculations 
and FOCUS [4] and Belle [2J] measurements. The intent 
was pedagogical, and we showed the plots at seminars 
and conferences (25| . 

Unfortunately, the error band in that effort was im- 
pressionistic, not rigorous. With the prospect of yet- 
more-precise results based on CLEO-c's full accumula- 
tion of 818 pb^^ [5], we now present a version that treats 
the error band as rigorously as possible. We also prepare 
plots for D — *■ 'kIv and B — > irlv. 

As before we shall base the plots for D decays on 
Ref. . The final result of this analysis consists of the 
BK parameters {F, a, [3) and the 3x3 error matrix. The 
full statistical error matrix is contained in a detailed, un- 
published description of a BK-based analysis of i? — > nlv 
form factors [2^. The best fit, statistical errors, and sys- 
tematic errors are tabulated in Table [H The statistical 
correlation matrices pij = CTfj/(o'^iO'|j)^^^ are tabulated 
in Table ini The correlations among systematic errors are 
discussed below. 

Propagating (correlated) fluctuations in F, a, and (3 
to the form factors, one flnds relative squared-errors 
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These errors are plotted as a function of in Fig. [T] as 
solid curves. The relative statistical errors are smallest 
for such that 



(13) 
(14) 



It is illustrative to take a"^^ and a|n^ from Tables HI and HIl 
and solve Eqs. (fT3| and (fT4|) for g^. We call these values 

and g| and tabulate them, as well as. qj^^^^^ and q^^y., 
in Table IIIII As one can see from Fig. [T| and Table [ml 
the statistical error is smallest between g^j^ g p-j and gi^axi 
as expected. One may view this outcome as a check on 
the fitting procedures. 

One can reverse this strategy to determine the corre- 
lation between the systematic errors of F and the slope 
parameters. In the error budget of Ref. @ the largest 
systematic effect comes from discretization errors. These 
should be smallest around q^^ q and g^ax for /+ and /o, 
respectively, because those correspond to the smallest |p| 
yielding the respective matrix elements. This yields 



syst 
PPa 



-0.198 {D 



Ppfi 



K), 



-0.329 {D 
-0.533 {D 



(15) 
(16) 



and the dashed curves in Fig. [TJ 

It is customary to combine statistical and systematic 
uncertainties by adding the two cr^ (matrices). Carrying 
out this procedure leads to the curves and bands in Fig. [21 
The error bands seem to contradict the conventional wis- 
dom that the lattice-QCD uncertainties are smallest near 
9max- This is not entirely the case for the relative error, as 
seen in Fig.[Tl As increases, the relative errors decrease 
until hitting a minimum somewhere between q^^ ^ and 

9maxj is reasonable. The form factors rise faster than 
the relative errors drop, leading to the increasing absolute 
error seen in Fig. [21 These features are not an artifact of 
the BK parametrization, as we shall see below with the 
B -kIv form factor /+. 

Figure[2lis the first main result of this article. It shows 
the form factors /+ and /o for D Klv and D irlv. 
The lattice-QCD results are shown as curves (red for /+, 
blue for /o) with two errors bands, one statistical (orange 
for /+, gray for /o), the other systematic and statistical 
combined (yellow for /-i-, light blue for /o). Experimental 



measurements for 



are overlaid as points 



with error bars. It may require careful scrutiny to see 
which experiment is which, but a glance reveals how well 
the points and curves agree. The agreement is good for 
D ttIv and very good for D Klv. 

For the z expansion the propagation of errors is even 
simpler. Focusing on /+, one has from Eq. ([5]) 



N 



^2 ^k+l 

O 1.1 ^ 



fi 



Z^k=0 "-k^ 



(17) 



where the indices on correspond to those on the se- 
ries coefficients. The coefficients and error matrix for the 



TABLE I: Best-fit values of BK parameters with statistical 
and systematic errors, successively, in parentheses [^. [tI [2^. 



Decay 



F 



D^Klv 0.73(3)(7) 0.50(4)(7) 1.31(7)(13) 
D-^ttIp 0.64(3)(6) 0.44(4)(7) 1.41(6)(7) 



TABLE IL Statistical error correlation matrices pij 
(^fj /{'^'ii'^'jjy^^ of the BK parameters |26j |. 
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-0.597 


1.000 
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FIG. 1: Relative errors vs . Solid (dashed) curves show 
the fitted statistical (estimated systematic) error for /+ (red 
curves) and /o (blue curves). Vertical lines show 9^1,0,0) 
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TABLE III: Useful quantities for generating and assessing Figs. [l][2l and|3] 



Decay 


H* 






2 

9(1,0,0) 


2 

9m ax 


-2 

9(1,0,0) 


~2 

9m ax 


-2 


9l 






(MeV) 


(MeV) 


(GeV^) 


(GeV^) 










D Klv 


Dt 


2112 


704 


1.10 


1.88 


0.25 


0.42 


0.47 


0.38 


D ~* nlf 


D* 


2008 


518 


1.57 


3.00 


0.39 


0.74 


0.53 


0.52 


B Tvlf 


B* 


5325 


518 


22.4 


26.4 


0.79 


0.93 







iV = 3 fit (to = 0.65g,^a^) are tabulated in Tablc[lVl The 
z-series fit was carried out after assigning g^-dependent 
systematic uncertainties, so Table HVl refers to the com- 
bined statistical and systematic errors of this analysis. 

This information, combined with the outer function 0+ 
0, is used to produce Fig. [31 the second main result 




FIG. 2: Form factors /+ (and /o) for semileptonic D decays, 
from lattice QCD 0, [2g, expressed as a red (blue) curve with 
an orange (gray) statistical error band and a yellow (light 
blue) combined error band. Error bands take correlations 
into account. Measurements of /+ are from Belle (green di- 
amonds) '2^, BaBar (magenta squares) (SJ, and CLEO-c 
(maroon triangles) ,2^, i2£] . The vertical line shows 



TABLE IV: Best-fit values at and correlation matrix pki of 
the 3-term z expansion of /+ for B nil/, with statistical 
and systematic errors combined 



Fit: 


0.0216(27) 


-0.0378(191) 


-0.113(27) 


P 


ao 


0.1 


a2 


fflo 


1.000 


0.640 


0.475 


ai 


0.640 


1.000 


0.964 


a2 


0.474 


0.964 


1.000 



of this paper. Now the curve and error band conform 
with preconceptions, for several reasons. First, q^^ ^ is 

close to Qn-jaxi rather than in the middle of the kinematic 
range. Second, the chiral extrapolation in Ref. is less 
aggressive than that in Ref. [a], leading to a larger but 
more realistic error at 5^ = 0. The most striking aspect 
is that even though the absolute error in is increasing 
for > 0.8, the band remains narrow. The band simply 
conveys the point-to-point correlations better. 

This paper shows in detail how to compare semilep- 
tonic form factors from lattice-QCD and from experi- 
ments. For illustration we use the BK parametrization 
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lattice QCD [Fennilab/MILC, 08 1 1 .3640 [hep-lat]] 

■ experiment [BaBar, hep-ex/06 12020] 
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FIG. 3: Form factor /+ for B —> nlv expressed as a curve 
(red) from the best fit with a total error band (yellow) from 
taking correlations in the fit parameters into account over- 
laid with measurements of \Vub\f+/{3.38 x 10~^ 
(magenta squares) [30|]. The vertical line shows g, 



from BaBar 

2 

max ■ 
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for D Klv and D ttIu, and the z expansion for 
B ttIv. Clearly, the idea is more general. For ex- 
ample, an interesting prospect relevant to semileptonic 
form factors is to inject 3-momenta smaller that P(i,o,o) 
using "twisted" boundary conditions [3l|, [H, [s^l . That 
strategy should improve the accuracy of parameters in 
the chiral extrapolation and, hence, the BK, BZ, or z 
fits. The output of any fit could still be exhibited as 
outlined here, although one should bear in mind that su- 
perior visualization of a fitting procedure does not repair 
any shortcomings of the fit itself. 
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